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The possible role of growing static order in the dynamical slowing down towards the glass tran¬ 
sition has recently attracted considerable attention. On the basis of random first-order transition 
(RFOT) theory, a new method to measure the static correlation length of amorphous order, called 
“point-to-set (PTS)” length, has been proposed, and used to show that the dynamic length grows 
much faster than the static length. Here we study the nature of the PTS length, using a polydis¬ 
perse hard disk system, which is a model that is known to exhibit a growing hexatic order upon 
densification. We show that the PTS correlation length is decoupled from the steeper increase of 
the correlation length of hexatic order, while closely mirroring the decay length of two-body den¬ 
sity correlations. Our results thus provide a clear example that other forms of order can play an 
important role in the slowing down of the dynamics, casting a serious doubt on the order agnostic 
nature of the PTS length and its relevance to slow dynamics, provided that a polydisperse hard disk 
system is a typical glass former. 


Introduction 

When we supercool a liquid while avoiding crystalliza¬ 
tion, dynamics becomes heterogeneous mm and slows 
down significantly towards the glass transition, below 
which a system becomes a non-ergodic state. Now there 
is a consensus that this slowing down accompanies the 
growth of dynamical correlation length [3]. Several dif¬ 
ferent physical scenarios have been proposed, yet the ori¬ 
gin is still a matter of serious debate: while some sce¬ 
narios describe the glass transition as a purely kinetic 
phenomenon [4], others posit a growing static order [5] 
or a loss of configurational entropy [ 6 ] behind dynamical 
slowing down. Among this last category, we will focus 
here on two distinct approaches. The first one is ran¬ 
dom first-order transition (RFOT) theory [iHS], which 
is based on a finite dimensional extension of mean-field 
models with an exponentially large number of metastable 
states. The second approach, recently proposed by some 
of us [TOl [11] , ascribes the growth of the dynamical corre¬ 
lation length with the corresponding growth of the static 
correlation length. Here we focus on these two scenar¬ 
ios based on static order and consider which is more 
relevant to the origin of glassy slow dynamics, using a 
simple model glassformer, two-dimensional (2D) polydis¬ 
perse hard disks [I1[I3|. 

In RFOT theory, metastable states are thought to have 
amorphous order, whose correlation length diverges to¬ 
wards the ideal glass transition point. It was recently 
suggested that the so-called point-to-set (PTS) length, 
which is the correlation length of amorphous order, can 
be extracted by pinning a finite fraction of particles and 
studying the dependence of the overlap function on the 
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pinning particle concentration. According to the RFOT 
theory, amorphous order develops in any glass-forming 
liquids and this method is thought to be able to pick up 
the static correlation length whatever the order is, i.e., 
the method is claimed to be order agnostic m- Thus, the 
use of pinning fields has been considered to be a promis¬ 
ing new direction in the study of the glass transition. 
Within the RFOT theory, it was shown that freezing the 
positions of a finite concentration of particles shifts the 
ideal glass transition to higher temperatures, potentially 
granting access to the glass state in equilibrium [T5H22]. 
Moreover, the average distance between pinned parti¬ 
cles at the liquid-to-glass transition represents a direct 
measure of the PTS correlation length. PTS correla¬ 
tion lengths aim at measuring hidden static length scales 
by looking at the extent of the perturbation induced by 
frozen particles on the rest of the liquid. It is intuitively 
defined as the average distance between pinned particles 
that forces the system to stay in an amorphous config¬ 
uration with a vanishing configurational entropy. The 
reasons behind the popularity of PTS correlation lengths 
in the study of the glass transition are at least twofold: I) 
they are expected to provide an “order agnostic” method 
to measure static correlations [SHE!]; 2) it is theoreti¬ 
cally established that no divergence of the relaxation time 
of a glass at finite temperature can occur without the con¬ 
comitant divergence of the static correlation length [25] . 

On the other hand, it was recently noted nniiiiiiisi 
that, for moderately polydisperse hard disks, an increase 
in the area fraction of particles 0 , hexatic (or 6 -fold bond 
orientational) order grows and its correlation length is 
supposed to diverge, obeying the Ising-like power law, to¬ 
wards the ideal glass transition point ^o, where the struc¬ 
tural relaxation time diverges following the Vogel- 
Fulcher-Tammann law. We have also confirmed that the 
dynamical correlation length ^4 is proportional to the 
hexatic correlation length ^6 and furthermore there is al¬ 
most a one-to-one correspondence between the degree of 
hexatic order and the slowness of dynamics. These re- 
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suits suggest an intimate link between static order and 
dynamics: the dynamical slowing down is accompanied 
by an increase in both size and lifetime of hexatic or¬ 
dered regions. We also found that 3D poly disperse hard 
and Lennard-Jones spheres exhibits essentially the same 
behaviour m- These results suggested that the dynam¬ 
ical slowing down is a consequence of the growing acti¬ 
vation energy associated with the Ising-type power-law 
growth of the correlation length of critical-like fluctua¬ 
tions of static order towards the ideal glass transition 
point [iQiin]. Recently, a theory for the occurrence 
of such criticality in disordered systems with topologi¬ 
cally ordered cluster of particles was also proposed by 
Langer [26]. 

The role of local order on the dynamics of glassy sys¬ 
tems remains controversial at least for two reasons. The 
first problem is that the local order that one needs to 
measure is system-dependent, and up to now the rele¬ 
vance of bond orientational order was demonstrated only 
for polydisperse particle systems and a spin liquid, and 
not for bidisperse systems mm- The second problem 
is conceptual: are static correlations really responsible 
for the dynamical slowing down? The PTS correlation 
length is often described as a remedy to both problems, 
since it should be able to detect static correlations with¬ 
out a detailed knowledge of the local order involved in 
these correlations. In studies of binary mixtures of hard 
spheres, the PTS correlation length was shown to grow 
only modestly in the regime accessible to computer sim¬ 
ulation Haul], differently from the dynamical correla¬ 
tion length which grows much more rapidly. These re¬ 
sults suggested that no link exists between a single static 
length scale and the dynamical slowing down (the only 
exception would be close to a possible ideal glass transi¬ 
tion temperature) [5l|T4]. On the other hand, measures 
of the PTS correlation length [28] have shown that it cor¬ 
relates well with the average dynamics of the system, and 
with dynamic heterogeneities [29], The PTS correlation 
length is in agreement with measures of the density of 
plastic modes [30], providing support for the idea of a 
fundamental length scale controlling the dynamics of the 
supercooled liquids. 

Unlike previous studies, in this work we measure the 
PTS correlation length in a system for which a growing 
local order was previously found, i.e., polydisperse hard 
disks. This will allow us to compare the growth of PTS 
correlation lengths with the correlation length of bond 
orientational order. We will consider several pinning 
strategies (random pinning, uniform pinning and cavity 
pinning) (Fig. [T -C), and then look for the underlying 
structural features that are captured by the PTS corre¬ 
lation length. In principle, each different pinning geom¬ 
etry probes a different lengthscale [31]. The point-to-set 
static lengthscale ^pxs was first introduced in the spher¬ 
ical cavity geometry [32] . Both random and uniform pin¬ 
ning are expected to express the same static lengthscale, 
here called and the RFOT theory predicts a temper¬ 
ature scaling relation between ^pxs and which in its 


simplest form is written as ^k{T) ^ ^pt5(F)^/^ [I?]. If 
the PTS correlation method is indeed order agnostic, it 
should be able to pick up the correlation length of hex¬ 
atic order in a polydisperse hard disk system, provided 
that it is a typical glass former. Thus, it is a main inter¬ 
est of this work to reveal whether the PTS length is the 
same as the hexatic correlation length. This question is 
of crucial importance to reveal the origin of slow glassy 
dynamics. 


I. RESULTS: UNPINNED CASE 

The system studied is composed of N = 10000 polydis¬ 
perse hard disks with disk-size polydispersity A = 11% 
(see Methods). We start by considering the case without 
an external pinning field, c = 0. We focus on the follow¬ 
ing number densities from p = 0.92 to p = 0.97, which 
correspond to area fractions ranging from (j) = 0.73 to 
0 = 0.77. As described in Supplementary Information^ 
we extract the correlation length for bond-orientationally 
ordered regions by fitting the exponential decay of the 
peaks of the correlation function ge{r)/g{r) (see Fig. SI). 
We use an exponential function instead of a 2D Ornstein- 
Zernike function to avoid a priori assumptions on the 
origin of the growth of the correlation length. The corre¬ 
lation length ^6 is plotted in Fig. IDE, together with other 
length scales which we will derive later. The two-body 
correlation function ^2 is obtained by fitting with an ex¬ 
ponential law the decay of g{r) — 1. The results of this 
fit are also summarized in Fig. [T]E. These results confirm 
that, for polydisperse glass forming systems, the growth 
of many-body correlations associated with bond orienta¬ 
tional order is much faster than the growth of two-body 
correlations [nun Eg. 

As was shown in Refs. mm there is a link between 
dynamic heterogeneities and regions of high hexatic or¬ 
der in the fluid. In Fig. [T]D we show the results of Event- 
Driven molecular dynamics simulations in the isoconfig- 
urational ensemble (see Methods), where 200 trajectories 
are started from the same initial configuration but with 
different initial velocities. The degree of structural order 
is investigated by taking the average hexatic field over 
these N = 200 trajectories (also called isoconfigurational 
average) at a time t = Tq,/10 , where is the struc¬ 
tural relaxation time measured through the intermediate 
scattering function. The dynamics is instead investigated 
through the isoconfigurational average of the relative dis¬ 
placement |Ri| after t = Tq,, which is also approximately 
the time at which the heterogeneities are maximum (as 
measured by the four-point susceptibility [3]). The rela¬ 
tive displacement R^ is defined as the displacement be¬ 
tween time t = 0 and t = u(^a) ~ ^(O): of particle i 

with respect to its M neighbours, R^ = ^ Vj. 

This operation is introduced since coherent translational 
motion of a group of particles does not contribute to the 
stress relaxation. All disks in the configuration of Fig. [T]d 
are then grouped into sets of high and low mobility/order 
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FIG. 1: Types of pinning fields and the growth of static and dynamic correlation lengths, (top row) Different 
pinning strategies, with pinned particles coloured in black: random pinning (a), uniform pinning (b) and cavity pinning (c). d, 
Snapshot of a configuration at p = 0.97 in which the disks are coloured according to the following criteria: White^ low mobility 
and high order; Black, high mobility and low order; Cyan, low mobility and low order; Magenta high mobility and high order, 
e, Correlation lengths as a function of density p: red squares for bond orientational order ^q, blue triangles for the two-body 
correlations ^ 2 , black circles for random pinning, green diamonds for for uniform pinning, black filled circles for ^pts- 
The inset shows the scaling of the structural relaxation time Ta with the hexatic correlation length (points), and the fit with 
the relation Ta = tq exp(T)^/(fo) (dashed line). 


depending whether their mobility/order is higher or lower 
than the 50th percentile. We can then identify four differ¬ 
ent sets of particles: low mobility and high order (white); 
high mobility and low order (black); low mobility and low 
order (cyan); high mobility and high order (magenta). 
Our results show that 76% of particles are in the first two 
sets (38% in each), demonstrating a high degree of corre¬ 
lation between structural ordered regions and immobile 
regions (or, vice versa, between disordered regions and 
mobile regions). Moreover, the remaining two sets (each 
accounting for the 12% of particles) are located at the in¬ 
terface between mobile and immobile extended regions. 
In Fig. magenta disks are located on the surface and 
in between black clusters, while cyan disks are located on 
the surface and in between white clusters. In other words 
disks which are next to an low/high mobility region, will 
also have low/high mobility. Here we note that the em¬ 
bedded fractal nature of order parameter fluctuations is 
characteristic of critical fluctuations. 


II. RESULTS: PINNING 

Having characterized the static properties of the un¬ 
perturbed system, we now introduce the pinning field. 
First, simulations are fully equilibrated with cluster- 
moves algorithms (see Methods). A representation of 
the pinning fields is given in Fig. [^- C. In the random 
pinning geometry, Np particles are chosen randomly and 
pinned, i.e. their position is kept fixed during the course 
of the simulations. For each density we introduce pinning 
fields with concentrations c = 0.01,0.06,0.10,0.15,0.20, 
and for each concentration we average over 9 different 
realizations of the fields. In this scheme the distance be¬ 
tween pinned particles is defined only as an average over 
a broad distribution, as both clusters of pinned particles 
and extended regions without pinned particles are likely 
produced. For this reason, random pinning is expected 
to be more sensible to finite size effects, as was observed 
in Ref. [19], where it is noted that random pinning can 
smear out the Kauzmann transition in very small sys¬ 
tems. The pinning geometry can also have strong effects 
on the dynamics [lUISlJIllIll]. In order to limit the 
fluctuations in the distance between pinned particles, we 
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also adopt a uniform pinning geometry, where a simple 
cubic lattice is overlaid to the equilibrated configuration, 
and the closest particle to each lattice point is pinned. 
Particles are pinned at the following average distances: 
a = 2.5,2.75,3,4.25,6,8,10, and each distance is aver¬ 
aged over 7 realizations of the field. Finally, in the cavity 
pinning geometry, all particles outside a cavity of ra¬ 
dius R are pinned. Since the static length scales cur¬ 
rently accessible to simulations are expected to be small, 
well within lOa, simulations with cavity pinning involve 
a small number of particles, thus requiring extensive av¬ 
erage over different realizations of the field (here 100 sim¬ 
ulations for each cavity diameter). 

Since the pinning field is applied to equilibrium config¬ 
urations, the static properties should be unchanged with 
respect to the c = 0 case. We check this by computing 
both positional and hexatic order for different concentra¬ 
tions c. All results are consistent with the c = 0 case and 
the standard deviation between simulations at different c 
is represented with the error bars for in Fig.[2 Corre¬ 
lation lengths are extracted from all pinning geometries, 
following the procedure outlined in Supplementary Infor¬ 
mation. In all cases the physical idea is to detect the char¬ 
acteristic length (average distance between pinned par¬ 
ticles in the random and uniform pinning geometries, or 
the size of the cavity in the cavity geometry), which pro¬ 
duces a high localization of the mobile particles, as mea¬ 
sured by overlap functions. We plot the random pinning 
correlation length random, the uniform pinning corre¬ 
lation length uniform, and the point-to-set lengthscale 
from cavity pinning ^pts in Fig- a® (see Supplementary 
Information on the details of its estimation). We see that 
the growth of the PTS lengthscale, irrespective of the 
pinning strategy, is significantly slower than the growth 
of bond orientational correlation length ^6, while being 
comparable to the growth of pair correlations in the sys¬ 
tem ^2 • We also confirm that estimating the PTS length 
from uniform and random pinning through the relation 
fpTs{T) ^ Ci'(F) still produces a much weaker growth 
than that of ^6- The inset of Fig. shows the scaling 
between the ^6 length scale and the relaxation time r^, 
= ^0 exp(I)^/(^o), where D is a measure of the fragility 
of the system m- This scaling also supports a direct 
connection between the growth of structural correlation 
and slow dynamics. 

All static lengthscales can also be obtained by con¬ 
sidering coarse-grained variables, dividing the simulation 
box into smaller boxes of side length I = 0.3cr, ensur¬ 
ing that each box can be occupied at most by one disk 
at any time during the simulation. The problem is then 
mapped on a set of discrete variables, defined as nf = 1 
when the center of a disk is in box i in configuration cf, 
and nf = 0 otherwise. The average occupancy for a par¬ 
ticular realization of the pinning field is defined as nl, 
where the overline denotes an average over thermal fluc¬ 
tuations for a fixed realization of the pinning field. Figure 

shows the field fii for three different concentrations of 
the pinned particles for the case of uniform pinning at 


p = 0.95. Going from the low to high concentration (left 
to right in the figure) the amount of localization in the 
field progressively increases. Localization is associated 
with regions with high occupancy probability, such as 
the purple dots that can be seen for average distances 
a = 3(7 and a = 2.5(7. At a = 3(7, these high occupancy 
regions are localized in particular regions of the simula¬ 
tion box, while at a = 2.5(7 they fill the box rather uni¬ 
formly. This suggests that the PTS correlation length, 
which characterizes the crossover of the localization tran¬ 
sition, is between 2.5(7 < ^pts ^ 3 (7, which is consis¬ 
tent with the measurement based on the spatial correla¬ 
tion of the overlap function, ^pts(p = 0.95) « 2 .6(7 (see 
Fig.[T]E). Moreover we confirm that all static lengthscales 
obtained through the coarse grained representation are 
fully consistent with those obtained in Fig. [T (see Sup¬ 
plementary Information). The similar scaling behaviour 
between pinning lengthscale fx and the pair correlation 
length ^2 can also be understood in terms of the nl field. 
Each pinned particle generates an oscillatory perturba¬ 
tion of the rq field, which originates from the two-body 
static correlations between the pinned particle and the 
mobile particles in the liquid. So no localization tran¬ 
sition can occur if the average distance between pinned 
particles is bigger than the range of two-body correla¬ 
tions, ^ 2 - The localization transition thus requires pinned 
particles to be within ^ 2 , below which the number of par¬ 
ticle arrangements drastically decreased, and the config¬ 
urational entropy vanishes. The extent of these regions 
with high localization of particles is exactly what is being 
measured by the PTS correlation function. 

The results thus show that the growth of the pinning 
correlation length is similar to the growth of two-body 
correlations. On the other hand, the growth of the bond 
orientational order correlation length is much faster, and 
clearly decoupled from the pinning correlation length. 
We also checked that the same is true for coarse-grained 
quantities. This strongly indicates that the PTS correla¬ 
tion length is not order agnostic^ but targets the growth 
of a particular order in the system, that is the size of the 
regions where particles are localized due to the pinning 
field. The growth of these regions follows the growth of 
two-body correlations: particles are localized due to the 
perturbation that pinned particles introduce to the nl 
field, and the length scale of this perturbation is given 
by two-body correlations. In other words, at least in the 
density range considered here, the localization transition 
due to point pinning requires that the average distance 
between pinned particles is smaller than the two-body 
correlation length. 


III. RESULTS: LOCALIZATION 

Next we consider whether the localization transition 
is linked with the underlying hexatic ordering. A visual 
inspection of the configuration at a = 3 in Fig. shows 
that the localized regions tend to form rather compact 
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FIG. 2: Development of rii field with a decrease in a for uniform pinning. The particle density is p = 0.95, and 
the average distance between particles a = 6 (left), a = 3 (middle) and a = 2.5 (right). The colour code is mapped to the 
occupancy probability according to the colour bar. Distances are in units of the average diameter (a). 


domains. These domains are found by plotting the prob¬ 
ability distribution function for the occupancy field as 
shown in Fig. ||4. The figure plots P{n) for different val¬ 
ues of a, the average distance between pinned particles, 
at density p = 0.97. Without pinned particles (a = oc) 
the distribution shows one Gaussian peak centred around 
n = where I = 0.3cr is the coarse-graining length. For 
finite values of a, the distribution progressively broadens, 
and, for a < 3, covers almost all the n range, n G [0; 1]. 
A similar behaviour is also seen for the distribution func¬ 
tion of the coarse grained hexatic field -06, as shown in 
Fig. [^, where the distribution of is plotted for the 
same state points of the top panel. A visual inspection of 
the underlying configurations indeed confirms that local¬ 
ized particles appear in regions where the average hexatic 
field is strong, thus revealing an important structural fea¬ 
ture of the localization transition. Localization appears 
in regions of high hexatic order, where the average dis¬ 
tance between pinned particles is within the two-body 
correlation length, ^ 2 - 

Now we examine whether the localization observed in 
our simulations is a true first order transition or just a 
crossover. Within RFOT theory, the lower critical di¬ 
mension is two, so a true thermodynamic transition is 
not expected in our simulations, even though the ques¬ 
tion still remains open. Both distribution functions in 
Fig. and B do not show signs of bimodality. We next 
examine the global fluctuations of the overlap function, 
P{q) = {6{q — qap)) (see Fig. for p = 0.97). The av¬ 
erage value of q changes continuously as a function of a, 
and no sign of bimodality is present. For small values of 
a and high q the distributions deviate significantly from 
the Gaussian shape and display a heavy tail. This is not 
associated with the glass-transition localization but to 
a different localization process which happens when the 


distance between two pinned particles is small enough to 
block the passage of fluid particles. This is analogous 
to the mechanism responsible for the Lorentz-gas transi¬ 
tion [36] . We thus find no evidence of a first-order local¬ 
ization transition. Our results of course do not exclude 
the possibility that such transition could occur at higher 
densities. The broad distribution observed in Fig. HH 
and B are in fact compatible with the presence of two 
population in a mixed state. The inset of Fig. plots 
the probability distribution function is plotted for a = 3a 
and for different values of p: the different curves cross at 
an isosbestic point, which is a common feature of bimodal 
distribution functions, meaning that the overall popula¬ 
tion becomes more localized as p is increased (the pres¬ 
ence of the isosbestic point is independent of the value 
of a). But simulations at higher volume fractions require 
considerablely larger system sizes (due to the rapid in¬ 
crease of the hexatic order correlation length, .fe), and 
are outside the scope of the present investigation. 


IV. DISCUSSION AND CONCLUSIONS 

We have extracted several static length scales from 
systems of polydisperse hard disks with polydispersity 
A = 0.11, in the range p G [0.92; 0.97]. The results con¬ 
firmed that the length scale associated with bond orien¬ 
tational order grows more rapidly than the length of pair 
correlations [niiiaES]. The use of pinning fields enabled 
the calculation of the PTS correlation length, showing 
that it grows only moderately with increasing supercool¬ 
ing, a result which is in agreement with measures of the 
PTS length in binary mixtures [23] . For polydisperse sys¬ 
tems, the PTS correlation length is not coupled to that of 
bond orientational order, which is directly linked to the 
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FIG. 3: Dependence of the probability distribution functions of various fields on the average distance between 
pinned particles a for uniform pinning, a, Probability distribution function of the occupancy field ^ at p = 0.95 and 
different values of a. The inset shows P(n) for constant a = Scr and different values of p. b, Probability distribution function 
for the average hexatic field magnitude, iV^el, at p = 0.95 and for different values of a. c, Probability distribution of the overlap, 
P(^), for p = 0.95 and for different values of a. 


dynamical correlation length unmi]: the growth of the 
former is considerably slower than the latter. For differ¬ 
ent glass forming systems, this suggests that also other 
forms of order originating from many-body interactions 
could go undetected by PTS measures. 

The PTS length captures a localization transition that 
occurs in presence of pinned particles. This localization 
transition originates when the occupancy field, n, has 
extended regions of high probability due to neighbouring 
pinned particles. The perturbation that a single pinned 
particle produces in the n field is due to pair correlations. 
In absence of strong nonlinear effects, the length scale of 
the localized regions extends no further than two-body 
correlations. This is the case in the density interval ac¬ 
cessible to our simulations, where pinned particles need 
to be placed closer than the pair correlation length in 
order to produce localized regions in the fluid. A second 
requirement that our results suggest is that pinned par¬ 
ticles should be in positions compatible with high local 
hexatic order. 

The localization transition that occurs with increasing 
concentration of pinned particles happens continuously 
in the density range we could access in equilibrium. The 
results do not rule out the possibility that strong non- 
linearities will produce a localization transition that ex¬ 
tends beyond pair correlations for higher (but yet un¬ 


reachable) densities. 

To summarize, the PTS length measured by particle 
pinning simply reflects pair correlation and fails in de¬ 
tecting the correlation of bond orientational order (more 
precisely, hexatic order), which intrinsically originates 
from many-body interactions. Although the PTS length 
is decoupled from the dynamical correlation length, the 
hexatic order correlation is strongly coupled to it. This 
implies that slow dynamics in our system is controlled by 
the development of hexatic ordering, and not by transla¬ 
tional order detected by the PTS correlation. Although 
the generality of this conclusion needs to be checked care¬ 
fully, our study suggests that the PTS length is not order 
agnostic and the growth of the PTS length is not respon¬ 
sible for glassy slow dynamics at least for our system. 
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Supplementary Information 


Simulations 

We study two-dimensional polydisperse hard disks 
with Monte Carlo simulations. The diameter a of the 
disks is extracted from a Gaussian distribution, and the 
polydispersity is def ined as the s tandard deviation of the 
distribution, A = (cr)- In the present work 

we fix A = 11% for which no transition to an hexatic 
phase is observed. The unit of length is set by the aver¬ 
age disk diameter (a). 

All simulations are run at fixed densities p = N/V^ 
with N = 10000, with the event-chain algorithm mm, 
which allows for a fast equilibration even at very high 
densities. After the equilibration run, we activate the 
pinning field and switch to Metropolis dynamics with 
swap moves between randomly selected pairs of non- 
pinned particles. 

The connection between static and dynamic length- 
scale shown in Fig. ID was obtained with Event Driven 
simulations [3] in the isoconfigurational ensemble [4], 
where 200 trajectories are started from an equilibrated 
configuration at p = 0.97 but with a different assignment 
of initial velocities. 


Estimation of and ^2 

Here we explain how to estimate the hexatic order cor¬ 
relation length ^6 and the translational correlation length 
^2- Figure]^ shows snapshots of configurations at differ¬ 
ent densities, where each particle is coloured according to 
the phase of the local hexatic order, arg('06). The figures 
show that the hexatic order increases with p. 

Bond orientational order in two dimensional polydis¬ 
perse hard disks is expressed by the hexatic order param¬ 
eter 

tAe = — W 

k 

where rij are the neighbours of particle j, and Ojk is the 
angle that the bond between particle j and k makes with 
a reference axis. In the above definition, neighbouring 
particles are defined as particles sharing an edge in the 
radical Voronoi diagram obtained from the particle’s po¬ 
sitions and sizes. To measure the extent of hexatic order 
we define the following correlation function 

geir) = (V’6(»’)V’6(0)) 

and extract the correlation length from the spatial decay 
of g6{r)/g{r). 

Figure 1^ shows the hexatic order correlation function 
gQ{r)/g{r) for different densities. The correlation length 


is obtained by fitting the exponential decay of the peaks 
of the correlation function. We use an exponential func¬ 
tion instead of a 2D Ornstein-Zernike function to avoid 
a priori assumptions on the origin of the growth of the 
correlation length. 

Positional order is expressed by the pair correlation 
function, p(r), and its length scale measured by the de¬ 
cay of g{r) — 1. The two-body correlation function ^2 is 
obtained by fitting with an exponential law the decay of 

g{r) -1 - 

Estimation of random and ,uni form 

The static lenghscale was measured in both the ran¬ 
dom and uniform pinning geometries following the pro¬ 
cedure defined in Ref. [5]. First a microscopic overlap 
function is defined as 

Wmn{0,t) = 0(a - fn{t) “ rm(0)|), 

where 0 is the Heaviside function, and a = 0.3cr en¬ 
sures single occupancy of cells with side a. The overlap 
function is then averaged both over different equilibrium 
configurations and for different realization of the pinning 
field as 

Qc{t) = n _ Q\]\f ^ ^ 

m,n^V 


p = 0.94 p = 0.95 


p = 0.96 p = 0.97 

FIG. 4: Snapshots of configurations of polydisperse 
hard disks at different densities. Particles are coloured 
according to the phase of the hexatic order, with similar 
colours sharing the same orientation of the hexatic order. The 
size of the correlated regions increases with density. 
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FIG. 5: Decay of the hexatic order correlation function 

gQ(r)/g(r) for different densities. Continuous lines are ex¬ 
ponential laws obtained by fitting the peaks of the correlation 
function. 


where V is the set of pinned particles with concentra¬ 
tion c = Np/V, and the brackets stand both for thermal 
and disorder averages. The function Qc{t) is one at t = 0 
and then decays to an equilibrium value that is a mea¬ 
sure of the overlap between the initial configuration and 
the subspace of configurations that are compatible with 
the pinning field. For low concentrations c, the function 
Qc{t) is expected to decay rapidly, while for high con¬ 
centrations c it should decay to a high value. This cor¬ 
responds to two states: the fluid state with low overlap, 
and the glass state with high overlap (even in the ideal 
glass state the overlap is not Qc(oo) = 1 because of ther¬ 
mal fluctuations). By measuring Qc(oc) — Qc=o(^) 
a function of the concentration c, the point-to-set (PTS) 
correlation length is defined as ^ps = (cpsp)^^^, where 
c* is the concentration that locates the transition. For 
small system sizes, it was shown that, at low enough 
temperatures, this transition appears to be a first-order 
phase transition [6]. For big system sizes, the barrier 
which separates the two states become too high to be 
simulated, so one defines the PTS correlation length as 
the crossover length between small and large overlap. We 
follow Ref. [5] and set Qcpsi^) — Qc=o(oo) — 0.4. Care 
has to be taken in ensuring that the simulations are prop¬ 
erly equilibrated, which can be checked by the decay of 
the self part of the overlap function: 

QcW = (1 Wnn{0,t))- 

For polydisperse hard disks with 11% polydispersity, this 
limits the simulations to the range p < 0.97. 

Figure plots the average overlap function (5c(^) at 
the highest density considered p = 0.97. The function 
Qc{t) decays to a plateau, whose height depends on the 
concentration of pinned particles. For low concentra¬ 
tions, the overlap between any two configurations in the 
system is low, and Qc(oc) — Qc=o{^) decays to a low 



FIG. 6: Average of the overlap function Qdf) for dif¬ 
ferent concentrations. Here p = 0.97. 



FIG. 7: Overlap as a function of the average distance 
between the particles for different densities. The PTS 

correlation length is defined as the length at which Qcps (oo) — 
Qc=o(oc) = 0.4. For the lower densities this value is obtained 
by fitting points with < 10 by a polynomial and 

extrapolating the result. Given the rapid decay of the overlap, 
the results do not depend sensibly on the fitting procedure. 


value; for high concentrations of pinned particles, the 
system is constrained in a region of phase space with 
high overlap (and overall low configurational entropy), 
and Qc(co) — Qc=o(co) decays to a high value. The PTS 
correlation length is defined as the length that charac¬ 
terizes this crossover. To extract it, we follow Ref. [5] 
and plot the value of Qc(oc) — Qc=o(oc) versus the aver¬ 
age distance between pinned particles, (cp)“^/^, as shown 
in Fig. We fix the crossover value to be around 
Qc^{oo) - Qc=o{oo) 0.4, and define ^ps = [ckpY^!^^ 
We have checked that the results do not depend sensibly 
on the exact choice of the crossover value, as the overlap 
decay rapidly at the crossover. 

The same results are confirmed with coarse-grained 
variables. In this approach the simulation box is divided 
into smaller boxes of side I = 0.3cr, ensuring that each 
box can be occupied at most by one disk at any time 
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FIG. 8: Average occupancy rn at density p — 0.95. Here 
we show the n field for different concentrations of the pinning 
field, c = 0,0.01,0.10,0.20. 


during the simulation. The problem is then mapped on 
a set of discrete variables, defined as nf = 1 when the 
center of a disk is in box i in configuration a, and nf=0 
otherwise. The average occupancy for a particular real¬ 
ization of the pinning field is defined as nl, where the 
overline denotes an average over thermal fluctuations for 
a fixed realization of the pinning field. 

In Fig. we plot ni at density p = 0.95 for concen¬ 
trations c = 0,0.01,0.10,0.20. For zero concentration, 
the nl is structureless, with a uniform value of nl Pp^ 
where I is the side length of the boxes. When a pinning 
field is acting on the system, nl measures the effects of 
the perturbation. In the figure, pinned particles are rep¬ 
resented by black circles. Each pinned particle generates 
an oscillatory perturbation of the nf field, which origi¬ 
nates from the two-body static correlations between the 
pinned particle and the mobile particles in the fluid. This 
can be seen by measuring the spatial correlation function 
of the rii field, 

92 ^ in - rj) = {Wnj)p - {rH)p{rij)p, 

where the (• • • )p is an average over different realizations 
of the pinning field. By fitting the decay of 92 ^ {r)/g{r) 
with an exponential function, we extract the correlation 
length ^ 2 ^. This length scale is equivalent to the two- 
body correlation length of the unperturbed system (with¬ 
out pinning field), ^2- 

The coarse grained PTS length is estimated from the 


following overlap function. 



where are boxes that do not contain pinned parti¬ 
cles, a and ^ are two arbitrary configurations, and the 
brackets denote both thermal and pinned field averages. 
We evaluate the PTS length of the coarse-grained field 
as described for the non-coarse grained variables, with 
(5^(oo) — Q^£q{oo) = 0.25. We have extracted the cor¬ 
responding lengthscale, verified that it matches 

exactly the ^pp^ lengthscale. 

The use of coarse graining is not limited to the oc¬ 
cupancy variable, nf. We can define a coarse grained 
hexatic field as 

a 

where Nt is the total number of configurations and 'ipQ j 
is the hexatic order for particle j in configuration a, and 
Sij = 1 if particle j is in box i, and zero otherwise. 
Pinned particles will produce perturbations also in the 
field. By taking the spatial correlation function, 
{pQ,i (^)'06,2(O)): we obtain the coarse grained correla¬ 
tion length, which, as expected, is the same length 
scale as measured by the bond orientational order in the 
unperturbed case, ^6- 


Estimation of ^pts 

The PTS lengthscale ^pts is measured in the cavity 
geometry following the procedure defined in Ref. [7]. A 
coarse grained description of the overlap function is intro¬ 
duced, in which space is divided N in cubic boxes of side 
a = 0.3(7. One then defines the occupancy number ^ 
which is 1 if box i in configuration and a is occupied and 
0 otherwise. The coarse-grained overlap is then defined 
as both thermal and disorder averages of the following 
quantity as 


1 ^ 

quit) = + 1))- 


The decay of ^ oc) = q{R) is then fitted with a 
compressed exponential function as 


where rj is the compressed exponent, and A and a are 
fitting parameters. 

An example of cavity pinning field is given in Fig. 

In the cavity pinning geometry, the overlap is strongly 
affected by the location of the cavity in the starting con- 
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From the simulations we extract the PTS length (see 
Methods) by studying the overlap in 7x7 boxes of side 
a = 0 . 3(7 located at the center of the cavity. The value 
of the measured overlap as a function of cavity size R is 
plotted in Fig. 

The data in Fig. [To] were fitted with a compressed ex¬ 
ponential: q{R) = A exp with A = 0.5, fol¬ 

lowing the procedure described in Ref. [7]. 


Intermediate scattering function 


FIG. 9: Snapshot of a configuration at p = 0.97 with a 
cavity of radius R = 9(7. Pinned particles are coloured in 
black, while particles inside the cavity are coloured according 
to their index. 

figuration: cavities inside a region of high hexatic order 
lead to higher overlaps, while cavities in disordered re¬ 
gions are characterized by smaller overlaps. A large num¬ 
ber of realizations of the pinning field are thus required: 
in our simulations we average over ten different initial 
configurations equilibrated at density p, and for each con¬ 
figurations we sample ten cavities in random positions in 
space, for a total of one hundred simulations for each 
density p = 0.92,0.4,0.96 and cavity size R/a G [2,9], 
with AR = ±1(7. 



FIG. 10: Overlap (symbols) for cavities of size R/a G 
[2,9], with AR = ±1(7. The results are shown for densities 
p = 0.92, 0.4,0.96,0.97. The lines are best fits to a compressed 
exponential (see text). 


Figure plots the intermediate scattering function, 
Fs{t), for different p values, showing an increase of the re¬ 
laxation time by more than two orders of magnitude, and 
a relaxation that acquires a two-step stretched exponen¬ 
tial character. Observing the dynamics at p > 0.97 re¬ 
quires considerably larger system sizes, as the correlation 
length of the hexatic order parameter becomes compara¬ 
ble to the simulated system size {N = 10000 disks) (see 
Fig. !§. The increase of the relaxation time (shown 
in the inset) follows the usual Vogel-Fulcher-Tammann 
(VFT) law Tq, = Aex.p{—Dp/{p — po)), where best fits 
give D = 0.42 and po = 1.018. 



FIG. 11: Density-dependence of self-intermediate scat¬ 
tering function, Fs{t). It is computed at the wavelength 
corresponding to the first peak in the structure factor for den¬ 
sities in the range p G [0.92,0.97]. The inset shows the struc¬ 
tural relaxation time, Ta computed as Fs{ra) = 0.1 (sym¬ 
bols) and the VFT fit (line), = Aexp(— F)p/(p —po)), with 
D = 0.42 and po = 1.018. 
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